Anderson-localization versus derealization of interacting fermions in one dimension 
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Using the density matrix renormalization group algorithm, 
we investigate the lattice model for spinless fermions in one 
dimension in the presence of a strong interaction and disorder. 
The phase sensitivity of the ground state energy is determined 
with high accuracy for systems up to a size of 60 lattice con- 
stants. This quantity is found to be log-normally distributed. 
The fluctuations grow algebraically with system size with a 
universal exponent of ~ 2/3 in the localized region of the 
phase diagram. Surprisingly, we find, for an attractive inter- 
action, a delocalized phase of finite extension. The boundary 
of this delocalized phase is determined. 

PACS numbers: 71.30.+h, 72.15.Rn 



The influence of electron-electron interaction on Ander- 
son localization has attracted a lot of interest for several 
years. Many recent studies were motivated by the exper- 
imental observation of persistent currents in mesoscopic 
rings Motivated by an early suggestion Q that the 
interaction between the electrons may give a significant 
contribution to the average persistent current, this phe- 
nomenon in the presence of both interaction and disorder 
has been investigated by various methods |J-§|. Never- 
theless, the magnitude of the effect is still not well un- 
derstood. 

In one dimension, interacting systems in the absence of 
disorder [| ll|, as well as for disordered systems in the 
absence of interactions [jl2j are well studied. However, 
a clear understanding of the interplay between interac- 
tion and disorder has not yet been obtained. In this 
Letter, we present novel results of a detailed study of 
a simple interacting-fermion model with disorder. We 
determine the ground state phase sensitivity with high 
accuracy for a wide range of parameters and system sizes 
up to 60 lattice constants. Our main results are (i) a 
universal behavior of the rms-value of the logarithmic 
phase sensitivity, which grows with system size, M, pro- 
portional to Af 2 / 3 in the localized region, and (ii) the 
zero-temperature phase diagram, which shows, for an at- 
tractive interaction a delocalized phase of finite exten- 
sion. 

The numerical results are obtained with the density 
matrix renormalization group algorithm (DMRG) [ p"3| ], 
which allows calculation of ground state properties of 
disordered, interacting fermion systems with an accu- 
racy which is comparable to exact diagonalization, but 



for much larger systems |L4|,[L5||. In our implementation 
of the DMRG we perform 5 finite lattice sweeps keeping 
up to 750 states per block. 

We consider a chain of spinless fermions with nearest- 
neighbor interaction and disorder, 
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and twisted boundary conditions, c = e^CM- The 
length of the chain is denoted by M, and the particle 
number is N. For simplicity, we will set t = 1 in some of 
the formulas below. 

The ground state energy E(<j>) depends on the phase <j>. 
The energy difference between periodic and anti-periodic 
boundary conditions, AE = (~) N [E(0) - E(w)], the 
persistent current, I(4>) ~ —E'(<j>), and the charge stiff- 
ness, D ~ E"((f> = 0), are a measure of the phase sensi- 
tivity of the system. In the clean limit, i.e. e„ = for all 
n, the ground state energy can be determined from the 
Bethe Ansatz |^,[l7) . At half filling, the phase sensitivity 
in the limit of large systems (M — * oo) is given by ]I6| ] 
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where Em is the ground state energy of the M- 
site system, Soo is the energy density in the ther- 
modynamic limit, and v is the Fermi velocity, v = 
■nt sin(2f7)/(7r — 2rj). The interaction parameter is K = 
7r/4?7, where r] parameterizes the interaction according to 
V = -2tcos(2ri). Thus MAE = nvK/2. 

For the noninteracting system, the phase sensitivity in 
the presence of a single defect (eo ^ 0) can be determined 
easily |L8) . In the presence of both electron-electron in- 
teraction and defect, it is more difficult to calculate AE. 
However, it is known that the ground state and the low- 
lying excitations can be described within the framework 
of a Luttinger liquid pT| , p9| . Combining a first-order per- 
turbative calculation with the scaling equations of Kane 
and Fisher |20| , we find in the case of a weak impurity 
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where Mo is a short distance cut off which is - for the half 
filled band - of the order of the lattice spacing. Using 
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a duality relation between a weak impurity and a weak 
link |20| , we obtain for a strong impurity 
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since the transmission through a strong defect is pro- 
portional to £ 2 / e o- Figure [j] shows MAE as a function 
of interaction for system sizes M = 60 and several de- 
fect strengths. The points are numerical results from 
the DMRG, the full lines are analytical results from the 
equations above. The cut-off parameter Mq was fitted in 
order to obtain agreement between numerical and ana- 
lytical results, giving Mq ~ 2. An attractive interaction 
makes the barrier more transparent, while a repulsive in- 
teraction increases the defect strength. Deviations from 
the analytical results are found when the "strong" impu- 
rity becomes so weak that a first order expansion is no 
longer appropriate. In addition, we find further devia- 
tions near V = ±2, where the Luttinger liquid becomes 
unstable. At V = —2, there is an instability with respect 
to phase separation (v — > 0, K — > oo). At V = +2 there 
is an instability to formation of a charge-density wave, 
since at this point 4/cj7-backscattering processes become 
relevant. 

In the presence of a weak random potential, but V = 0, 
we find, generalizing the single impurity result, 
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where we introduce disorder by taking the {e n } uniformly 
random distributed over the interval [— W/2, W/2]. At 
half filling, kp = 7r/2, the sum ^2 n=1 e„(— 1)™ can be con- 
sidered as a onedimensional random walk with M steps. 
Recalling that a random walk leads to a Gaussian dis- 
tribution of distances, we obtain for the average phase 
sensitivity and the fluctuations 
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The brackets (• • ■) denote the impurity average. It is 
apparent in Eq. (^) that perturbation theory breaks 
down for arbitrarily weak disorder, if the system is large 
enough: as is well known, in 1-D even weak disorder leads 
to localization with a localization length which is propor- 
tional to t 2 /W 2 . In large systems, MAE is drastically 
reduced due to disorder, however, numerical results (us- 
ing exact diagonalization methods, which is straightfor- 
ward as long as V — 0) indicate that it remains positive 
for all realizations of the disorder [ |T5[ , in agreement with 
a theorem by Leggett ^lj . For large systems we find an 
exponential decay of the average phase sensitivity. In the 
localized regime, i.e. M > £, the logarithm of MAE has 



approximately a normal distribution jig] . From our nu- 
merical data, where we averaged over 10 4 realizations of 
the disorder potential, and considered systems of up to 
10 3 sites, we find for the average logarithmic phase sen- 
sitivity and its variance a 2 in the limit of large systems 
(M^oo): 
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with £ = 114i 2 /W 2 . In order to check the universality 
of the exponent, we calculated the phase sensitivity for 
strong disorder up to W — 15t, and for different fillings. 
We always found the exponent 2/3 in the localized region. 

The interaction changes some of the results described 
above drastically. Applying the Kane-Fisher scaling to 
Eq. (ph we obtain 
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since the strength of each defect is renormalized. The 
average phase sensitivity is then given by 
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Again, a repulsive interaction tends to enhance the effec- 
tive strength of the defects, and an attractive interaction 
reduces it. Especially, for K > 3/2, i.e. V < —1, the 
strength of each defect vanishes so fast that disorder be- 
comes an irrelevant perturbation: there is no localization 
p2|,p3[ . We discuss the localized phase, V > — 1, first. 
Assuming that only one relevant length scale exists, i.e. 
the localization length £, one concludes from ( p"l| ) that 
£ oc W 2 /( 2K ~ 3 ^ for weak disorder. This is verified in Fig. 
^ where we plot the logarithmic phase sensitivity as a 
function of the scaled systems size. In the case of the 
largest systems considered (M = 60), we averaged over 
several hundred realizations, whereas for short systems 
[M < 20) we used ensembles of more than 10 3 real- 
izations. With good accuracy, points corresponding to 
different strengths of disorder lie on the same curve, i.e. 
the localization length is indeed the only relevant scale, 
even for M 3> £, where perturbation theory breaks down. 
The average phase sensitivity, shown in Fig. ||, is for large 
systems approximately given by (V = 1.2) 

(MMAE)) =-M/£+l, (13) 

with the localization length £ sa 28W^ 2/( - 3 ^ 2K \ 



2 



The rms- value, cr\ n (MAE)i shown in Fig. |j| is for small 
systems proportional to M < ~ 3 ~ 2K ^ 2 , see Eq. (|l2|). (Note 
that <J\ n (MAE) and cmae are directly related to each 
other, provided <7mae "C (MAE).) For large systems we 
again find the fluctuations to be proportional to M 2 / 3 , 
as in the noninteracting limit. Explicitly, we found from 
our numerical data (V = 1.2, i.e. K 0.71) 
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In Fig. | we plot (\n(MAE)) as a function of interac- 
tion and for several system sizes (here W = 1). For com- 
parison, we included the phase sensitivity in the absence 
of disorder. Between V ~ —1.6 and w —1.1, the phase 
sensitivity remains almost unreduced, even for large sys- 
tems. We believe that this region corresponds to the de- 
localized phase predicted earlier 22 23]. This assertion 
is confirmed by an apparent divergence of the localiza- 
tion length when approaching the phase boundary from 
the localized side [g4|. Nevertheless, the phase sensitiv- 
ity remains smaller than in the clean system since the 
parameters v and K scale downwards due to the random 
potential p3| . 

The fluctuations of the logarithmic phase sensitivity 
provide another, more accurate method for determining 
the extension of the delocalized phase. Selected data are 
shown in Fig. || for M = 10 and 30. As discussed in 
connection with Eq. (|l^), a decreasing variance (with in- 
creasing M, compare O with +) implies that the disorder 
scales to smaller values, hence the system is delocalized, 
while the variance increases for a localized ground state. 
Using this property as the criterion, we obtain, with con- 
siderable numerical effort, the phase diagram shown in 
Fig. [| These results are based on system sizes between 
30 and 50 sites. Clearly, by this method, we can only 
give a rough estimate of the phase boundary, and it is 
possible that we somewhat overestimate the size of the 
delocalized region. 

In summary, using the DMRG algorithm, we have ob- 
tained high accuracy results for the ground state energy 
for a model of interacting fermions with disorder. In 
the weak disorder limit, we verified quantitatively sev- 
eral predictions on disordered Luttinger liquids. In the 
localized region, we determined the localization length 
and the distribution of the phase sensitivity. The lat- 
ter is nearly log-normally distributed, with a universal 
size dependence of the fluctuations proportional to M 2 / 3 . 
We confirmed the existence of a delocalized region in the 
phase diagram. As far as we know, we are the first to 
give a quantitative estimate of the size (as a function of 
disorder and interaction) of this region. 
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FIG. [T|: Phase sensitivity of the ground state energy in 
the presence of a single defect as a function of interaction 
for several defect strengths eo- The diamonds and crosses 
are numerical results (system size M = 60). The straight 
lines are analytical results as described in the text. 
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FIG. Average logarithmic phase sensitivity as a func- 
tion of the scaled system size, for V = 1.2 and disorder 
W = 1,2,3. 
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FIG. |3|: Rms- value of ln(MA_E) as a function of scaled FIG. |5|: Rms-value of the logarithmic phase sensitiv- 
system size (again V — 1.2). The full line is the analytic ity versus interaction, for M = 10 and 30; W = 1. 
result according to Eq. (12), which explains the low-M 
behavior. For large systems, cr\ n (MAE) is proportional to 
M 2 / 3 (dashed line) as in the noninteracting case. 
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FIG. U Average logarithmic phase sensitivity as a func- 
tion of interaction, for system sizes ranging from 10 to 
40; W = 1. For comparison, we included the result in 
the clean limit (W = 0, dashed-dotted line). 
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FIG. § Phase diagram. The symbol O (+) denotes the 
region where the variance of the logarithmic phase sen- 
sitivity decreases (increases) as a function of the system 
size. We considered up to 50 sites. The O-region corre- 
sponds to a delocalized ground state. 
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